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Abstract 

We determine the general form of the potential of the problem of 
motion of a rigid body about a fixed point, which allows the angular 
velocity to remain permanently in a principal plane of inertia of the 
body. Explicit solution of the problem of motion is reduced to inversion 
of a single integral. A several-parameter generalization of the classical 
case due to Bobylev and Steklov is found. Special cases solvable in 
elliptic and ultraelliptic functions of time are discussed. 
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1 Introduction 

Integrable systems constitute a rare exception in mechanics. This is most 
clearly manifested in the field of rigid body dynamics, where all known 
integrable cases constitute few small tables, see e.g. p], [2] and [3]. Second 
comes particular solutions of equations of motion of a rigid body in various 
settings. Those are solutions valid only for certain particular sets of the 
initial position and angular velocity. In the classical problem of motion of 
a rigid body about a fixed point in a uniform gravity field there are eleven 
solutions of this type known after authors of the 19th and the 20th centuries. 
All of them are collected in table 1 below (in chronological order) : 
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Table 1: Known particular solvable cases of the classical problem. 

For a detailed account of those cases see [16] or [T7] . Some of them were 
generalized through the addition of a gyrostatic moment [16] and other 
potential and gyroscopic forces [E], [IB] . 

In the present article we aim at exploring the possibility of particular 
solutions of the Bobylev-Steklov type for the problem of motion of a rigid 
body about a fixed point in a field that generalizes the classical setting. We 
assume that the body is acted upon by certain potential forces, which admit 
symmetry axis fixed in space. The equations of motion for this problem can 
be written in the Euler-Poisson form (e.g. pQ): 
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= 73^— 
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dV 


= 7i^— 
012 


~ 72-5— 
071 



Bq + (A - C)pr 



71 + 973 - ^72 = 0, 72 + rji - P73 = 0, 73 + p-/ 2 - <37i = 0, (2) 

where A,B,C are the principal moments of inertia, p,q,r are the compo- 
nents of the angular velocity of the body and 71,72,73 are the components 
of the unit vector 7 along the axis of symmetry of the force field, all being 
referred to the principal axes of inertia at the fixed point. The potential V 
depends only on the Poisson variables 71 , 72 , 73 • In the classical problem of 
a heavy body V = 071 + b 72 + C73. 
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Equations (pQ) and ([2]) admit three general first integrals: 

I\ = -^^P 2 2^ r2 ^ ^ e ener ^ integral (3) 

1% = 7i + t| + 73 = 1) the geometric integral (4) 

^3 = Ap^i + -B(772 + Cr73, the areas integral (5) 

2 A new solvable case 

The aim of this paper is to look for potentials V which allow full solution of 
the system (Q3[2]) under the condition 

q = (6) 

on the zero level of the areas integral, i.e. 

h = (7) 

We assume that i / C. without restriction on the third moment of inertia 
B. This choice of the problem is motivated by the classical Bobylev-Steklov 
solution in the classical problem, which is characterized by the same condi- 
tion ([SD but with the potential V = a"fi and the additional restriction on 
the moments of inertia A = 2C. 

From ([6]) and the middle equation of (P) we get 

dV dV 

(A-C)pr- l3 — + 11 — = (8) 

Differentiating this equality and using equations ([I])-© we arrive at the 
expressions for the remaining angular velocities 



/ ^73 , dV_, 



—^—(-n— -73— ) 0) 



where V satisfies the linear partial differential equation 



3 



. A d 2 v dV. ,, 2 _ 2 , . d 2 v d 2 v 



9 7 f 5 7 | 572^73 071^72 

<9F <9F <9V" 

-04 - 2C) 7 273^ + 2(4 - C)^— - (24 - C) lll2 — = QflO) 
071 1772 073 

On the other hand, from the first and third equations of ([5]), in virtue of 
([6]), we have 

*73 _ 73 _ CJ3 

dryi 71 4 71 
This can be readily integrated to give 

73 = ^7i^ A , A = const. (12) 

This means the trace of the vertical unit vector on the unit sphere attached 
to the body lies on a cylindrical surface whose generators are parallel to the 
y— axis. Also, from (J3|) one can express 72 as 



72 = Vl-7i 2 -A 2 V (13) 

asd Thus, five of the six Euler-Poisson variables are expressed in terms of 
the last one: 71. The relation with time can be determined by separation of 
variables in the first Poisson equation. This finally gives 

dryi 



Vain) 



(14) 



where 



/ \ A li A ° 2 x2 9V dV. 

9M = c(A-c)x {1 -^- x ^ )(7i ^- 73 5t7 )o 

and ()o in the right hand side means the value of the expression in virtue of 
relations (fl~2l [13j) , so that g is a function of the single variable 71 . 
Summing up, we formulate the following 

Theorem 1 For an arbitrary rigid body moving about a fixed point while 
acted upon by forces with potential V satisfying U0\) . the Euler-Poisson 
equations (QHSf are solvable on the zero level of the areas integral under the 
condition q = 0. The solution is parametrized in terms of 71 by expressions 
(0 EH 1 13\) and relation with time is given by 
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3 The general form of the solution 



It is not hard to construct the general solution of the linear PDE (jlOp . which 
may be written in the form 

V = V 1 + V 2 , (15) 

where 



Vl _ (All+Cli) I H m±jii^ ^ (16) 



V 2 = (A 7l 2 + C7 3 2 )G(72) (17) 

This form involves two arbitrary functions F and G, which we assume well 
behaved on the poisson sphere (jl]), so that all subsequent operations on the 
potential V are justified. With this form of the potential we can rewrite the 
expressions for the Euler-Poisson variables in their final form parametrized 
by 7i 



(p,q,r) = (-\J^-^i /A ,0,\l^-li)V^ 



(71,72,73) = (7i,Vl"7i 2 -A 2 V, A 7l 7M ) (18) 



where 



F {-\^(A - C)*" 1 ) f A^+CX^ A m ( 2 + A 2 2C M) _ u] A/C du 
™ = : ^7771 1" / F( ; ^7777 ) 



A>fi + C\* 7 f /A J C( 7l 2 + A 2 7 r M )-n 

I 2C 

+G(Vl-7i 2 -A 2 7 1 " x ), (19) 
and 71 is determined in terms of time by inverting the integral 
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d 7i , x _ ^7f 



2C 



, »(7i)=-2^(l-7i 2 -AV)^ (20) 

V9(7i) G 

It may now be verified that the solution given by the expressions (|15H20|) 
satisfies the Euler-Poisson equations ([H [2]). Also, one can check that the 
total energy of the motion is in fact preserved and has the value 

h = -F{-\-c{A-C)c- v ) (21) 
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4 Generalization of the classical Bobylev-Steklov, 
Chaplygin and Goriatchev cases 

It is difficult to foresee the explicit form of the part (|16p of the potential with 
a given choice of the function F in the integrand. It may also be impractical 
or even impossible to evaluate the integral in (|16p for arbitrary A and C 
in closed form. In the classical Bobylev-Steklov the moments of inertia are 
subject to a single condition A = 2C. We now explore the form of the 
potential corresponding to certain simple forms of the function F under this 
same condition, to obtain a generalization of the classical Bobylev-Steklov 
case. It turns out that the sequences of functions {u n+5 }, {u~ n } lead to 
rational potentials. As an example we take the first three terms of each, so 
that 

B B B 

F{u) = V^l{A 1 + A 2 u + a 3 u 2 ) + — + -J + — ^ (22) 



u u 2 u 3 



gives the potential 



Vi = ai 7l + a 2 [87i( 7l 2 + 7 |) + ^] 

7i 

+a 3 [16 7l (8 7l 4 + 16 7l 2 7 l + 9 7 f) + ^ + f 3 \ 

A^^ ^-'l^ 1 71 (23) 

7s 7s 7 3 

di, hi being arbitrary parameters. If also we choose G( 72 ) = C\ + C2 7 | and 
use the geometric relation @, we may write the full potential (|15p as 

V = V x + Cl ( 7l 2 - 7 |) + c 27 2 (2 7l 2 + 7 |) (24) 



The solution of the equations of motion corresponding to this potential 
may be expressed by inserting it into Q and then substituting expressions 
(fl~2"1 [T3|) for 7 3 , 72 • One can obtain 71 as a function of time by inverting the 
hyperelliptic integral 



2 

1 



g(x) = — (1 — x — Ax) 

x{A 4 [2c 2 x 5 - 2(4a 2 - c 2 A)x 4 - 2(2a 2 A + c x + c 2 )x 3 - (ai - a 2 A 2 )x 2 ] 
-2(6i A 2 + 36 2 )x + 26 2 A} (26) 
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This integral becomes ultraelliptic when C2 = and elliptic if, moreover, 
0-2 = &2 = 0. In the last case the whole solution may be written as follows: 



V = ai7i + ci(7i - 72 ) + "4 



(27) 



73 





P 



(26 1 +a 1 ^ 2 7l + 2c 1 i/ 2 7 1 2 ),g = 0, 



r 



--(2&i + aii/ 2 7i + 2ci^ 2 7 2 ) 



(28) 



t = 




•71 



a/ (x 2 + ux — l)(26i + i/ 2 aix + 2z/ 2 cix 2 ) 



(29) 



Thus, the solution of the Euler-Poisson system (p]i2] ) with the three-terms 
potential (|2T|) can be expressed explicitly, for the body satisfying A = 2C 
and on the zero level of I3, in terms of elliptic functions of time. 

It should be noted here that the the solution (|2"7l does not depend on 
the second moment of inertia I?, which remains arbitrary. If in this solution 
we add the condition B = A, so that the the body has the Kovalevskaya 
configuration A = B = 2C, we obtain an intersection with three well-known 
integrable problems of rigid body dynamics: 

1. The full potential (j27j) with three arbitrary parameters a\,cx,bi was 
pointed out first by Goriachev in [10], but only under the condition 
A = B = 2C. 

2. When 61 = 0, we get a problem of motion of a rigid body by inertia 



in an ideal incompressible fluid (see e.g. [[19]) 



3. b\ = c\ = 0,Kovalevskaya's case of motion about a fixed point in the 
uniform gravity field. This case is general integrable, i.e. on arbitrary 
level of I3. 

Formulas (|29p give the explicit solutions of the first two cases under 
the additional restriction q = and the third under the two restrictions 
q = I3 = 0, in terms of elliptic functions of time. 
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5 Some closed-form and polynomial potentials 



Most problems of physical importance and all known completely integrable 
and solvable cases in rigid body dynamics are characterized exactly or ap- 
proximately by potentials of simple forms that are mainly polynomial or 
algebraic in the Poisson variables 71,72,73- It is quite interesting to iso- 
late possible solutions of polynomial or finite form of the linear PDE (|1U|) 
satisfied by V for arbitrary moments of inertia. The part V2 of (|15p is poly- 
nomial for any polynomial choice of (7(72). The part V\ as expressed by (|16|) 
depends only on 71 , 73 and hence satisfies the following particular version of 
equation (jTOj) : 



»7i $73 vjl dj3 



dV dV 
-(A - 2C) 13 — - (2A - C)^— = (30) 

Oil OJ 3 

We now assume a solution of this linear PDE as a combination of homoge- 
neous terms of the two variables 71 , 73 . After some trials to normalize the 
singularities of the resulting differential equations we arrive at the following 
form 

2 

V = Y^a vl v l9v {v), v = l + -^ (31) 
u Tl 
where the summation, with arbitrary constant coefficients a u , extends over 
the set of possible values of v that will be determined later. This leads to 
the differential equations for each of the functions g v 

„(!_„)_& + [(„_„+ __),, _ _ + 2| _| + (_ _ 2)g „ = (32) 

in which a = ^. This is a hyper geometric equation whose solutions for 
generic values of the parameters are 

9*1 = n~A-^^y~^v), (33) 



9u2 



2{a —1)2 

^ a-l + v (34) 
a — 1 



The second solution g U 2 contributes to (|3ip a term 

AM + llY'^ll + 7l) = A v (l - 72 2 )^ 2 - 1 (a 7l 2 + 7l) (35) 



S 



which is of the type f)lT[) . and can be better considered as included in that 
form. It thus remains to consider the decomposition of the potential part 
V\ ([16p in the form implied by (|15p and f|33|) . The typical term will now read 



Fl „ =7rF( -- 1 -- 5 _ r - + 2 ;„) (36) 



As we are looking for closed-form solutions we try to isolate cases when 
the hypergeometric series terminates, giving a polynomial expression. To 
that end, we note that the hypergeometric function satisfies the relation 



F(a, b; c; v ) = (1 - v) c - a - b F(c -a,c-b,c) 



and thus 



F(--,l- • L + 2-v) = (l-i,) 1+ 2(^T)F(2,l+ ; - ---+2;v) 

K 2' 2(o-l) 2 ' ; v ; v ' 2(a-l) 2 ' ; 

There are 3 obvious possible types of such closed- form expressions, in which 
the hypergeometric series truncates: 

1) a is a negative integer, while c is not a negative integer greater than 
—n. This gives only one non-trivial possibility v = 2, which leads to poten- 
tial of the type (fT7|) for a constant £(72). 

2) 6 is a negative integer —re, while c is not a negative integer greater 
than —n. This gives the sequence {u = 2 ^ a ~^ (n + 1), n = 0, 1, 2, ...}. The 
third index for this choice becomes a + 1 + (a — l)n, which is positive for all 
positive n. This choice invokes an infinite sequence of terms in the potential. 
The typical term in this case has the form 

F(2,-n ;« + ! + («- l)n; 1 + 4) ( 37 ) 



2(cm+a-l) v ' ' v ' ' 'v 2 

The numerator of this term is a homogeneous polynomial of degree 2n in 
the two variables 71 and 73. The power 2(cm + a — 1) in the denomenator 
is always positive. This sequence of potential terms is undefined at 73 = 
for all n. 

3) c — b is a negative integer — n, while a is not one of the n + 1 rational 
numbers 1 & = 0,1,..., n. This leads to the sequence {1/ = — 2(a — 

l)(n + 1), n = 0, 1, 2, ...} and to potentials of the type 

7 a . 7 2n F( _ n> _^ (n+1);1 _ n+ ^ ;1 + 23) (38) 

a a 7| 
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The product of the second and third factors in this expression is a polynomial 
of degree 2n. The behaviour of the first term differs depends on the two 
numbers a and n. 

The final form of the potential becomes 



V = (^7i 2 + C 73 2 )E^2 

n 

+7 f^ £ B n ^- 1/a) F(-n, -?^l( n + 1); 1 - n + 2±1- 1 + 2|) 

n 1 1 

E ^r7i 2 ^(2, -n ; a + 1 + (a - l)n; 1 + |) 
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